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The formalism of subdynamics is extended to the functional approach of quantum systems, and 
used for the Friedrichs model, in which diagonal singularities in states and observables are included. 
We compute in this approach the generalized eigenvectors and eigenvalues of the Liouvulle-Von 
Newmann operator, using an iterative scheme. As complex generalized eigenvalues are obtained, 
the decay rates of unstable modes are included in the spectral decomposition. 
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ON . I. INTRODUCTION 

For quantum system with continuous spectrum, the presence of resonances (small denominators) cause the failure 
of the usual perturbative methods for computing eigenvalues and eigenvectors of the time evolution generator. These 
■ difficulties have been considered as a manifestation of general limitations to computability for unstable dynamical 
m ! systems § §. 

In the work of the Brussels- Austin groups on Large Poincare systems, we find an algorithm to overcome the problem 
of small denominators, which are eliminated by a "time ordering rule". This is a rule for the regularization of the 
perturbation terms, which can be interpreted as a generalized boundary condition where terms corresponding to 
excitation processes are past-oriented, while terms corresponding to the de-excitation and emission of radiation are 
future-oriented Jl| ||. 

The construction provides a new type of spectral decomposition of the Hamiltonian operator. For the Friedrichs 
' model it was shown |l| that an appropriate mathematical framework for the time ordering construction is the theory 
C> of rigged Hilbert spaces of the Hardy class, formulated by A.Bohm and M.Gadella [Q || || 0. 

For mixed states, the construction algorithm is a generalization of perturbation theory based on the subdynamics 
decomposition of the Liouville-Von Newmann superoperator L (2) Through a non unitary transformation fi, the 
superoperator L is transformed into an operator Q = f2 _1 Lf2, which is block diagonal in the degrees of correlations. 
Qh, The perturbative method to obtain the spectral decomposition of the intermediate operator is regularized by 
(Jj ■ imposing the " time-ordering rules" , which in this case means that conserving or increasing of degrees of correlations 
is future oriented, while decreasing of correlations is past oriented j^]. This prescription explicitly incorporates 
irreversibility to microscopic theories. 

Usually, these perturbative algoritms are implemented in the so called "box normalization", in which the quantum 
system is assumed to be included in a box with periodic boundary conditions, the size of the box becoming infinite 
at some stage of the calculations. To perform this limit, it is necessary to consider volume dependent factors both for 
the diagonal components of the density operator and the observables. In this limit the recurrence time of the system 
jj] ■ is pushed to infinity. 

The diagonal singularities of operators for large quantum system was discovered by Van Howe [|l0| jll| [jl3| . 
At the same time, I. Prigogine and coworkers JlJ] [Jl5| |l6| |L7) emphasized the importance of states with diagonal 
singularity in non equilibrium statistical physics. 

Based in the pioneering work of I. E. Segal jl8|, I. Antoniou et al |l9| |^0| developed a formalism for quantum 
systems with continuous spectrum without the box normalization. 

The quantum states are junctionals over certain space of observables O. Mathematically this means that the space 
S of states is contained in O x . Physically it means that the only thing we can really observe and measure are the 
mean values of the observables O G O in states p G S cO x ((0) p — (p \ O)). This is the natural generalization of the 
usual trace of the product of the density operator by the observable (Tr(pO)) which is not well defined for systems 
with continuous spectrum. 

In this paper, we extend the theory of subdynamics to the case of quantum systems with diagonal singularities 
where, as stated in references |ll| [^0| ]2l|], the states are considered as functionals acting on the space of observables. 
The extended formalism is applied to the Friedrichs Model. 

In section II, we summarize the functional approach to quantum mechanics. The theory of subdynamics Q Js) Q is 
summarized in section III, and it is extended by us to the functional approach. In section IV, the extended formalism 
is applied to compute the generalized spectral decomposition and the time evolution for the Friedrichs Model. 



> 



i 



II. FUNCTIONAL APPROACH TO QUANTUM MECHANICS 



In the usual approach to quantum mechanics, a pure state is represented by a normalized vector of a Hilbert 
space TL. The observables of the system are represented by self adjoint operators acting on TL. The mean value of an 
observable represented by the operator O in a pure state represented by the vector \ty) is given by 

(O)* = (*|0|*> 

The time evolution of the state vector is given by the Schrodinger equation 

ij t \* t )=H\* t ), 

where the operator H : TL — >TL is the Hamiltonian operator of the system. Schrodinger equation has the solution 

|* t ) = e~ im |* > 

Mixed states have no well defined state vectors, but a probabilities p a (p a > 0, ^2 a p a = 1) of being in the pure 
states represented by normalized vectors |^a)- Therefore, the mean value of an observable O is given by 

(O) = 5> Q <vfgo|* a >. 

a 

The mixed state can be represented by the density operator p = 5Z a p a |^ r a)(^ r a|, having the following properties 

(O) = Tr(pO), Tr(p) = 5> a (* a |* Q ) = = 1 

a a 

As each vector |^ a ) evolves in time according to the Schrodinger equation, the time evolution of p is 

p t = e- iHt p e im 
and pt satisfies the Liouville-Von Newmann equation 

ij t Pt = Lp t (1) 

Lp=[H,p] (2) 

In a more general approach, the set of all possible observables of a quantum system is represented by an algebra O, 
while the possible states are represented by a set S of functionals acting on O (S C O x ). 

The mean value (0) p of the observable O in the state p is given by the value of the functional ponO, which we 
denote by (p\0) 

(0) P = (p\0) 

The last expression is antilinear in p and linear in O, i.e. 

(a lPl + a 2P2 \0) = aKp^O) + a* 2 {p 2 \0) (3) 

{p\a x Oi + a 2 2 ) = ai(p|Oi) + a 2 (p\0 2 ) (4) 

The algebra O is chosen to be an algebra of self adjoint operators on the vector space H, and as the mean value of 
the observables should be a real number, we impose the following condition on the states 

( P \o) = (p\oy, ifo = ot (5 ) 

The generalization of the concept of trace is 

Trp=(p\I) = l (6) 
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where I is the identity operator in the algebra O. 

For the time evolution in Heisenberg representation, the states are time independent, while the observables evolve 
in time according to 

t — e Ue 

The time evolution p t of the states in Schrodinger representation can be obtained from 

{0) t = ( Pt \0) = (po\e iHt Oe- iHt ). 

From the previous equation we obtain 

(j t Pt\0) = i( Po \e iHt [H,0] e- iHt ) = i(fit\[H,0]) 

Calling h'0 = [H, O] and omitting the observable O in the previous equation we obtain the generalized Liouville- 
Von Newmann equation 

-i(±p t \ = (p t p (7) 

To each superoperator M acting on S we can associate a corresponding adjoint superoperator acting on O and 
viceversa, with the formula 

(M P \0) = (p\M?0). (8) 
Using equation (||) with M = L, and the antilinearity property (j|) , in equation (]?]) , we obtain 

■ d TT 

which is formally the same as equation (|l|). However, as p is not an operator, equation (||) is not more valid. The 
extended Liouville-Von Newmann operator L acting on S is now given by 

<Lp\0) = {ppO) = (p\[H,0]). (9) 



III. THE FORMALISM OF SUBDYNAMICS 

Let us consider a linear space of states 5, and a linear operator L on S which is the generator of the time evolution 
of the states, i.e. 

i^-p t = hp t , p t eS. (10) 
at 

Let us assume that the operator L can be decomposed into 

L = L + Li (11) 

where Lo and Li are respectively called the " free" and " interaction" parts of L. It is also assumed that an interaction 
parameter A is included in Li to modulate the interaction. 

Starting with the projector Po on the invariant parts of the dynamics (LoPo = PoLo = 0), projectors P„ (n = 0, 1, ...) 
arc defined in such a way that they satisfy 

P„P„- = 5 nn ,P, l7 L P„ = P„L , P » = ^ P ™ ( L i)" Po = { J I I f n t 2 ( 12 ) 

n=0 3 

where I is the identity operator on S. The last equation means that the transition from F p to P m p is a process of 
m-th order in the interaction parameter. The operator P m es called the projection on the m-th degree of correlation. 
The main idea of subdynamics is to decompose the states through projectors II n (n = 0, 1, ...) satisfying 
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n »n„' = s nn ,Il n, Ln„ = n„L, n « = ^ = p "' (13 ) 

n=0 1_ * 

i.e. projectors n n commuting with L which reduce to the projectors P„ on the degrees of correlation when the 
parameter of the interaction tends to zero. 

Operators C n and B„, called creation and destruction of correlations are defined by 

c„ = q„c„p„ Q„n„ = c„p„n n 

B„ = P„B„Q„ n„Q„ = n„P„D n (14) 



where Q n = I — P n . From equations (|10|) and (|14j) it is obtained 

z^-(p„n„p) = e„p„n nP (15) 

at 

0„ = P„LP n + P„LC n P„ (16) 

n„ = (P„ + C„)(P„ + D„C„)- 1 (P n + B„) (17) 

[L , P m C„] = (P m C„ - P m )La(Pn + Cn) (18) 

[L ,B n P m ] - (P„ + D„)Li(P m - B„P m ). (19) 

The last two equations have the form [Lq , X] = Y, having the forward (backward) solutions 

±00 

X ± = i [ dt e-' Lot Ye ,Loi 



The following "time ordering rule" is chosen: +(-) sign is used for X = P m C„ with m > n (m < n), and for 
. = IDnPm with n > m (n < m). Therefore, the creation (destruction) of correlations, is future (past) oriented, i.e. 



±00 



±00 

-ihot 



= 1 / dte-^ ot (¥ n +B n )h 1 (¥ m ~B n P m )e lW m « (20) 
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The two previous equations can be solved iteratively to the required order in the interaction parameter, starting 
with the zero order solutions 

Cl°) = Bg» = 0. 

Once C„ and B„ are obtained, the intermediate superoperators 0„ = P„0„P„ can be computed using (Jl6|). The 
block diagonal super operator O = 0„ satisfies 

n 

l = nen~\ 



n = J2( F n + c «) . ( 21 ) 

n 

=Y / (Pn+ DnCn)- 1 (P„ + D n ) , 



and therefore it is isospectral with the Liouville-Von Newmann operator L. This property can be used to obtain the 
spectral decomposition of L in terms of the spectral decomposition of the superoperator O, with the same generalized 
eigenvalues. 

The formalism of subdynamics originally stated on the space of density operators can be easily translated to the 
functional approach of quantum mechanics described in the previous section. All the formulas of this section are still 
valid, but we should remember that the superoperators are in this case defined on the space S of functionals. 

It is operationally more convenient to rewrite all the previous equations for the corresponding adjoint superoperators, 
acting on the space of observables O. This is easily done by using the adjoint relation 
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(aMN^O) = (/>|a*N t M t O). 
For example, equations (|lj) for the projections on the degrees of correlations are replaced by 

^K' = 6 nn >K> 4^ = E P n = ^ ^(Ll)"FL = {;Jj^<~, (22) 

n 

where is the identify superoperator on 0. 

Equations (BO) for the creation and destruction of correlations transform into 



±00 



Ctpt, = -i / dte'^iPl + Ct)Ll(Ctpt l -P.L) e ^, m > n (23) 



< 







FlBl = -i / dte-^OP^-P^Dtjij^t +Dt) e a4*, m <n, (24) 



where now 



Ct=ptctQt, Dt = QtBtpt , Q^lt-P], (25) 



Equations ( J23| ) an (Q) can be written in a form which is more suitable for calculations. Let us consider the 
generalized left and right eigenvectors (a\ and \(3) of Lq, having degrees of correlation n a and np, i.e. 

(a|l4=o; a (a| ) Lj|/3)=^|/3), (a| = (a|pt Q , |/?) = P^ 

From equation ( p3| ) we obtain 

±00 

(a|Ct |/3) = -t y dte*^"^* (a|(pt +Ct)Lj (Cj, - Qt) |/3), n >n a 


If we use in the previous expression the identity 

±00 

( 



we obtain 



dt ~ ixt 



x±iO 



(ap n \l3) = ^— -(a\ (Pt + Ct ) Lj (C* - Qt ) |/J), n^n a (26) 

In the same way we obtain 

(a|Bt|/3) = L_(a|(Qt -B^L^pt + Bt)|/J), n^V (27) 

W/3 — W Q ± z(J > 

The intermediate operator 8' is 

©t e), = P^tpt + ptctLtpt (28) 

n 

The intermediate operator 0t is isospectral with LJ 

L f = (nt)-ietnt 

^ = E( p » + c «) ( 29 ) 



(Ot)-i=^( P t +Bt)(pt +C.tDt)- 1 
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For (PT 4- CJjB^) 1 we can use the following expansion 

oo 

(Pt + CtDt )-i = P t + £ ( C t pt )i (30) 

3=1 

The spectral decomposition of the intermediate operator 0t , i.e. a set of right (left) generalized eigenvectors \u na ), 
((u na \) satisfying 

(Una\u m p) = S nm S a p, 0^ = },Zna\Una)(Una\, P» = / ] \Una){Una\- (31) 

a a 

where a and (3 are discrete or continuous indexes. In the later case the sums in the previous expression should be 
replaced by integrals and the Kronecker by Dirac deltas. 

The spectral decomposition of LJ can be obtained from the spectral decomposition of the intermediate operator 
0t. From equations ( ^9|) and (|3l| ) it follows that 

= ^ ' ^nn |/na) (/to | j (32) 

where 

|/na) = (ft*)" 1 |«na) (/*a| = («na| « + - (33) 

The time evolution of a state functional, governed by the generalized Liouville-Von Newmann equation 

-il(p t \ = ( Pt \V, (34) 

is given by 

(pt\=Y, ei ' Z " a ' t (P°\f™)(fna\. (35) 



A. Perturbative solutions and A t approximation. 

If we replace the zero order approximation for n and lD>t (i.e. C^ '' = D^ '' = 0) in the right hand side of equations 
( p6| ) and (p^), we can obtain the first order approximations: 

(a|Ct«|/?) = =L_( a | pt L t(Qt ifl „ > (36 ) 

— w a ± it) < 

(a|B,t (1 >|/3) = ?__^( a | Qt L tpt |/3), ^< n „ (37) 

With Cn 1 ^ and it is possible to obtain the intermediate operators 0t up to second order, using equation ( |2S| ) 

0t( 2 ) = pt L tpt + pt CtC^L+pt (38) 
From equations (p9|), (|30|), ( |3^ ) and (^), we can compute fit and (f^) -1 up to first order 

nt(i) = £ ( pt +C t(i)) 

(n t)-KD = ^ (P t +D t(D) (39) 



The first order expressions for Ct and Ot given in equations (j^) and (|37| ) can be replaced in the right hand side 
of equations (Era) and (|27|) to obtain the next order approximation. In this way, through the computation of the 
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eigenvalues and eigenvectors of the intermediate operator 0t, it is possible to obtain the eigenvalues and eigenvectors 
of as a power expansion in the interaction parameter 



%na 


- z w 4 

not 1 


z {1) +z 

na 1 


\fna) 


= lf (0) ) 

1 J na J 




(/no | 


= (f (0) l 

\J na \ 


+ (f (1) l 

1 \J na \ 



F(2h 
na / 

f(2)| 



(40) 



Taking into account equations ([33|) and (|39| ) relating the eigenvectors of 0^ and we obtain 

= |5J2)+Di (1) l«J£) + + ••• 

(/na| = («<£| + (UJSIC^ + («&| + ■■■ (41) 

Replacing ( f40| ) and (|4l]) in ( |35j ) we obtain the following expression for the time evolution 

(pt\ = exp + z «« + z «« + •■•)* x 

na 

x(po\ [|5<2)(«W| + luWjC^ICJW + EtW|Sjg)(«(2| + |u«)(<>| + |u<2)(ii£>| +•••". 

If we omit first order contributions coming from the eigenvectors and third order contributions from the eigenvalues, 
the previous expression has the following approximated form 

( ft |^exp[i(^+,£)+a<] (Po\u^)(u^\. (42) 

na 

As we omitted first order terms in the eigenvectors, a necessary condition for equation (|42|) to be valid is 

A << 1, (43) 

where A is the interaction parameter. Moreover, as we considered the eigenvalues up to second order, a second 
condition involving the possible values of time is necessary for equation ( (42| ) to be valid 

X 3 t « 1, (44) 

which together with (^) gives 

t K \~ 2 . (45) 

In summary, the time evolution is given by the equation (^2|) if the interaction is small (A << 1) and the time is 
not too large (t<\~ 2 ). 

Conditions ([43) and (Eq) appear in the literature as the " X 2 t approximation". 



IV. FRIEDRICHS MODEL 

A. Observables, states and degrees of correlation 

The Hamiltonian of Friedrichs model is 

H = H o + V, 

OO 

H =m|l){l| +J du)U)\u))(w\, (46) 
o 

OO 

V = J duV u (\u){l\ + \l){w\). 



Where the vectors |1) and \ui) form a complete orthonormal set 
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(1|1) = 1, = (w|l) = 0, (w\uT) = 6 (w - u/) , /= |1)<1|+ / dw|a;)(w|. (47) 



o 



Let us consider the following definitions 

|1) = |1)(1|, |w) = \u)(u\, |1 W ) = |1)( W |, M) = | W )<1|, \u>«/) = \u>)(u>'\. (48) 
Any element O belonging to the space of observables O, can be written in terms of the operators defined in ( ft8| ) 

|0) = Oi|l)+ /" dwO w |w) + /" dwOi w |lw) + y <ia;O a , : L|wl) + y^y <&/CU'k<A (49) 

Notice that we explicitly included a diagonal singularity through the term J duiO^ui). 

We also consider the states as functionals acting on the operators. For this purpose it is convenient to define a set 
of functionals (1|, (u>\, (lw|, (wl| and (lolo'\ with the following properties 

(l\0) = O u (w\0) = O u , (lu\0) = lu , (ul\0) = O ul , (uu'\0)=O uu >. (50) 

or equivalently 



(51) 



In terms of these functionals, we assume that any element (p\ of the space of states S C O x can be written as 

(p\=P*i(M + I du J p* UJ {u J \+ f (L>p* lu) (loj\+ I dwp* 1 («l|+ jdu I du'plv, (uu'\ (52) 



(lll) = 


1, 


(1| W ) = 


(1| W 1) 


= (1|1 W ) 


= (l|u,a/)=0, 






(w\w') = 


5(w 




Mi) = 


(ul a;') 


= 1) = (uj 


) = 0, 




(lw\lw') = 


S(w 


-"'), 


(Mi) 


= (1 W | W ' 


= (lw|w'l) - 


(lwlw'w") - 


o, 


(wl\w'l) = 


5(u> 


-a/). 


( W i|i) 


= ( W 1| W ' 


= (wl|W) - 


(wl|a;V) = 


0, 


[u>u> \aa) = 


S(uj 


— a)<5(w 


'-a'), 


(wu/|l) 


= (ww' la) = ( 


u;u/|al) = (ww 



where 

Pl = Pl) pt=Pu, Pl u =Pul, Ptu>' = Pu'ui ( 53 ) 

P \+(<kjpl = l (54) 



Equations ( |53| ) are the conditions for p to be a positive functional, while (54) is a consequence of the total probability 
condition (p\I) = 1. In what follows (p\I) will be called the generalized trace of the state p (|J) = |1) + / da;|u;) is the 
identity operator in 0). 

By using the basis for states and observables, defined through and (0), we can also write 



lJ = / dw ( m — cj)|lCi»)(lo;| + / duo (uj — m)\cu 1)(uj 1| + J dui J dui' (u> — uj')\uj uj')(uj uj'\, (55) 
L f 1= y dwV r u [|a;l)-|lw)](l| + y dwV r u [|lw)-M)]( w l + y dw[-K„|l) + J du'V u >\u'u)](lu\ + 
+ MV.\D- du'V UJ ,\uu')](ul\ + J duj ydw'[K„|lw')-K'M)](ww'|, 

where Lj and Lj are the "free" and "interaction" parts of the Liouville-Von Newmann operator acting on O, i.e. 

LjO = [JJ 0l O] l\0 = [V,0] O eO (56) 
The diagonal and off-diagonal projectors, acting on O, are defined by 
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Pt = |l)(l|+ J cLj\uo)(uj\, 

ll = I duo\luo)(luo\ + I du>\wl)(wl\ + / duoduo'\uouo')(uouo'\ (57) 



The off-diagonal projector can be decomposed into 

QJ=p1+P|, F\ = J dL>\lu)(lu\ + J duj\ul)(ul\, Fl = J dw J duo'\uo u'\. (58) 

Pq, P{ and P\ are the projectors corresponding to degrees of correlation zero, one and two respectively, i.e. 
P\ (Li)°Pj=PjV0. 

pI(l\)°p\ ^pjpj = o, pJlJpJ^o. 

P$(lJ)°pJ=pJp£ = 0, Pjl4p|=0, Pj(Lt) 2 pt ^0. (59) 



B. Creation, Destruction and Intermediate Operators 



Using the equation (|36|) and (|37|) for the Friedrichs model we obtain for the creation and destruction operator up 
to second order are: 



4(1) 



da; 
duVu 
duo V u 
duj du> 



|1)( W 1| 



?71 — + iO CJ — TO + lO 

M)[(l|-M] IM[(1|-K 
— to + iO 



uo' — TO + i0 

I1 W )(1I 



to — w + iO 
\uo'uo){luo\ + — 



+ / dujduj' 



uo' — m — iO 



M)(l| 



TO — + lO TO — — iO 

I1)( W 1I ll)(M 



to — a;' + iO a; — to + iO 
uouo')(uol\ , 

wl)H 



+ 



UO — TO + iO TO — £J + i'O 
V U \1U)'){WU)'\ K;'|W1)(WW'| 



+ duo duo' Vu> 



\luo){uo\ 
m — + iO m — oj — iO 

\uouo'){uol\ |w'w)(1cj| 



to — a/ + iO a;' — to + iO 



a; — to + zO to — w + ?0 
Then, the intermediate operator Q^, up to second order, can be obtained using equation (|38| ) 

6j (2) =2^^|l)[(l|-(TO|], (m| = K 

/ dw[m-cj-/9]|lw)(lw|+ / du)[u — m-/3*]\wl)(ul\ + 
\luo')(luo\ \u'1)(uj1\ \oj'1)(1u\ 



duj dui Vui Vuj> 

duo duo' Vcj Vu> 
duoV 2 



TO — uo + iO TO — UO + «0 TO — cj' — iO 
|w'l)(wl| \luo'){uol\ \luo'){uol\ 



uo' — to + iO uo — to -f- iO uo' — m — iO 



Re(3 



duoV'P 



1 



Im[3 



UO — TO + zO 

duoduo'{uo — uo)\uo uo')(uo uo'\ + 

V u Vu" \uo"uo'){uouo'\ V u >V u »\ujuj")(u>u'\ 



-it ■ V 



duo duo duo" 



uo" — to + iO 



TO — zO 



(60) 



(61) 
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C. Generalized spectral decomposition and time evolution. 



From the explicit form of 0^ given up to second order in equations ( |6l| ) for Friedrichs model, the generalized 
eigenvectors and eigenvalues can be computed up to zero and second order respectively. The results are shown in the 
following table 









(^na | 


eS 


Zu> = 


\U U ) =J(w-m)|l) + \u) 
\ui) = \l) 


{u ul \ = (u\ 

(«i| = (1| - (ui = m\ 


el 


z lu = m — w — /3 

Zui = lu — m — [3* 


\Ul u ) = \lu>) 

\u u i) = \u>l) 


(Wlwl = (M 

(u„i\ = (oj1\ 


et 


Zulus' =U) - UJ' 


\u U u') = \^U>') 


Owl = 



(62) 



The generalized eigenvalues and eigenvectors given in the previous expression can be replaced in the equation d42j), 
to obtain the following time evolution 



(p t |l)-e- 2 ^*(p Q |l), 



(pt\w) = (p \u>) + 



1 - e 



-2nV„ t 



(po|l)«J(w - to), 



(pt\lco) 

(pM) 

(p t \LULu') = 



i{m—uj — (3) t 



i(uj—m—{3*) t 



(Po|lw), 
(PoM), 



(p \uJ) 



The first equation shows the decay of the discrete component (pt|l) of the state, with a rate Simultaneously, 
there is a growing term in the continuous distribution (p t \u>), with a sharp peak for the energy ui — to of the decaying 
mode. 



V. CONCLUSIONS 



We extended the formalism of subdynamics to the functional approach of quantum mechanics, in which the states 
are represented by functionals acting on the operators representing observables. 

The generalized spectral decomposition is obtained through an intermediate superoperator 0^ isospectral to the 
Liouville-Von Newmann superoperator (LJ = Q'0'O - The small denominators appearing in the perturbative 
expansions due to the continuous spectrum are regularized by the "ie-rule" (a time ordering prescription in which 
increasing (decreasing) of correlations is future (past) oriented). Due to this time ordering rule, 0t and therefore LJ 
may have complex eigenvalues. 

Considering eigenvalues up to second order and eigenvectors up to zero order, the time evolution is given by 



(Pt 



H ex P 

na 



where |un°2) and (u„a\ are generalized right and left eigenvector of ©J, computed up to zero order. For the previous 
expression to be valid, it is necessary that the interaction parameter be small and the time not too large, i.e. A« 1 
and t<\- 2 . 

When this procedure is applied to the Friedrichs model, we obtain 



(Pt\l) = e 
(p t |w) S (p |w) 



(Poll), 



(p \l)S(uj - to), 



(pM) 



(pt\vu>') = e 



i{m—uj — (3)t 



r-o i(u)—m—f3*)t 



(polM, 

(PoM), 



~ pi(u)-u)')* 



(po\wu)'). 



The first equation shows the decay of the discrete component (p*|l) of the state, with a rate Simultaneously, 
there is a growing term in the continuous distribution (p t |cj), with a sharp peak for the energy u — to of the decaying 
mode. 
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It is interesting to note that both decaying and growing terms are purely exponential. This may appear at first 
sight as a contradiction with the well known Zeno and Khalfin effects, which are deviations from exponential decays 
for small and big times. However, the previous expressions are not valid approximations for very big times. Moreover, 
Zeno effect implies 4t(pt\l)t=o- In our approximation, if we compute this derivative we obtain that it is of second 
order in the interaction parameter. As we neglected this order in the approximation, this result is not in contradiction 
with Zeno effect. 

In spite of the fact that the complex spectral decomposition can be obtained analytically for the Friedrichs model 
|, this paper shows that this approach is potentially suitable to deal with more complicated decaying 
processes,where it is impossible to obtain exact solutions 
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